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ON \ Abstract: Using a recently developed effective field theory for the interactions of nucleons 

at non-relativistic energies, we calculate the Coulomb corrections to proton-proton scattering. 
Including the dimension-eight derivative interaction in PDS regularization scheme, we obtain a 
modified Jackson-Blatt relation for the scattering lengths which is found to be phenomenologically 
satisfactory. The effective range is not modified by Coulomb effects to this order in the calculation. 
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\& , During the last year Kaplan, Savage and Wise have proposed an effective field theory 

for the interactions of nucleons at low energies These are characterized by scattering 
lengths which are much larger than the natural hadronic scale which in these systems 
is set by the pion mass. This made it initially difficult to construct consistent power 
counting rules which are necessary in any effective theory calculation to ensure that all 
contributions to a given order are included. By the introduction of a new regularization 
scheme called Power Divergence Subtraction (PDS) which is an extension of the more 
common MS scheme, these problems were solved. A very similar scheme was proposed at 
the same time by Gegelia[§|. Very recently, Mehen and Stewart |J have made this off-shell 



scheme (OS) more well-defined and shown that it is in fact equivalent to the PDS scheme. 

Within this framework several applications to problems involving the deuteron have 
now been madeQ^]. The first attempts to describe three-nucleon systems with two neu- 
trons have also been initiated 0. In this way the properties of the triton nucleus can 
be investigated. But also systems like 3 He with two protons are of obvious interest to 
consider in this new approach. For this to succeed, one must know how to include the 
repulsive Coulomb force between the two protons. As a first step in this direction we will 
here consider proton-proton scattering at very low energies. 

The effective Lagrangian for non-relativistic protons with mass M in the spin-singlet 
channel is 

A) = P f M^d + ^g^j P ~ ^-(P^2P)(p^2p) Jf (1) 

when we only include the lowest order contact interaction parameterized by the coupling 
constant Cq. It corresponds to the singular potential Co<5(r) which will affect interactions 
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only in the S-wave. In addition, we must also include the static Coulomb repulsion between 
the protons. It has the effective strength r\ = aM/2p where p is the CM momentum of 
protons, and thus also becomes strong at very low energies. This can be done using the 
well-established formalism based upon the exact solutions of the Schrodinger equation in 
the Coulomb potential Q. Denoting these by ip(r), the probability to find the protons at 
the same point in the initial or final state is reduced and given by ^(O)! 2 = C 2 where 

is the Sommerfeld factor 0. In a partial wave expansion of the full scattering amplitude^], 
the total phaseshift is ai + 5f where ag is the pure Coulombic shift and 8f is the shift due 
to the strong interactions. If we denote this latter one in the S- wave by 5p p , it is related 
to the corresponding strong amplitude Tsc(p) by the standard partial wave expression 

Plcoug -,) = --— (3) 

The superscript C is a reminder to the fact that this is not the strong phaseshift in the 
absence of Coulomb effects. Instead it is the strong phaseshift modified by the Coulomb 
potential. It is only the Coulomb repulsion between the protons in the initial and final 
states that have been removed at this stage. 

As has been known for a long time, cot 5p p in (|3j) does not have a regular effective 
range expansion Q. One finds instead 



P 



(cot 5g - i) + 2riH(rj)] = + K p 2 + . . . (4) 



n u 2 



where dp p and tq is respectively the S'-wave Coulomb-modified scattering length and effec- 
tive range. They arise after removing the part of the amplitude described by the complex 
function^] 

H (77) = il>{irj) + — ln(«7) (5) 

Lir\ 

which represents Coulomb effects at short distances. Its imaginary part cancels the unitary 
term ~ i in @. The real part defines the function h(r]) = Heip(ir]) — In 77 which is more 
suitable for phenomenological analysis! 11]. 

By this approach one were able to extract the two scattering parameters a long time 
ago, cipp = — 7.82fm and the effective range ro = 2.83 fm |ll]]. The scattering length is in 
magnitude much smaller than for pn and nn scattering in the spin singlet channel where 
one finds respectively the values a pn = — 23.75 fm and a nn = — 18.8fm. The difference 
between these two latter numbers is well within the range determined by the breaking 
of isospin invariance in strong interactions^]. For proton-proton scattering on the other 
hand it became clear that not all short-distance effects of the Coulomb potential were taken 
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out by the H- function in (||). By investigating the detailed behavior of the Schrodinger 
equation for the coupled problem of strong and Coulomb forces at small separations Jack- 
son and Blattfi^l succeeded in finding an expression for the scattering length a pp due to 
only the strong interactions, 
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0.33 



(6) 



With the above value for the effective range ro we then find a pp = — 17.0fm which is very 
close to the nn scattering length. The term —0.33 in the above formula is found to be 
only weakly dependent on the exact form of the strong potential p^j. 

From the effective Lagrangian ([!]) we can now calculate the proton-proton scattering 
amplitude. It will consist of the sum of all chains of proton bubbles separated by the 
vertex with value Co- Including the Coulomb potential corresponds to having static photon 
exchanges between all proton lines in these Feynman diagrams. The Coulomb interactions 
in the initial and final states just give the Sommerfeld factor (|2|). In perturbation theory 
each bubble will be replaced by an infinite sum of bubbles with zero, one, two and so on 
Coulomb interactions as shown in Fig.l. 





+ 





+ 



Figure 1: The Coulomb-dressed bubble is a sum of bubbles with static Coulomb exchanges. 

Since the scattering amplitude represented by the diagrams in Fig. 2 forms a geometric 
series, it thus takes the value 



Tsc{p) = C- 



C e 



"1-CoJoO) 



(7) 



Here Jo (p) is the value of one bubble which equals the amplitude for a pair of protons with 
energy E = p 2 /M to move from zero separation to zero separation. This is given exactly 
by the the Coulomb propagator Gc(E; r' = 0, r = 0). In momentum space it becomes 



Hp) = m 



d 3 k 2irr](k) 
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(2vr) 3 e 2 ^( fc ) - 1 p 2 - k 2 + ie 



(8) 



When this result for the scattering amplitude is now used in (^) and (Q), we see that both 
the phaseshift ctq and the Sommerfeld factor CqC 7 ?) cance l ou t- We are thus only left with 
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Figure 2: The proton-proton scattering amplitude is a sum of chains of Coulomb-dressed bubbles. 
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the evaluation of (||) which we do using the PDS regularization scheme[|lj in d = 3 — e 
dimensions and with a renormalization mass fi. It is ultraviolet divergent which shows up 
as an 1/e pole in the integral. This will be cancelled by counterterms which renormalize 
the coupling Co in (^) to Cq{u). As a result, the finite part of the dressed bubble @ is 
found to be 



j finite. n A/ 
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where Ce = 0.5772 • • • is Euler's constant. The last term here is the contribution from the 
special PDS pole in d = 2 dimensions. We see now that also the function H(rj) cancels 
out in (||). There is no contribution to the effective range vq to this order in the effective 
theory. Defining the /i-dependent strong scattering length 

1 4vr 



app(n) MCo(p) 



+ fi (10) 



as for pn and nn interactions^, we thus see from (||) and (^) that it can be expressed in 



terms of the measured scattering length 



' -pr + olM 



a pp{^) a pp 



(11) 



This relation has the exact same form as the Jackson-Blatt relation (|6|) where now the 
renormalization length 1/u enters instead of the effective range. But a pp {a) is not a 
physical quantity in the sense that it can be measured directly and will thus in general 
depend on the renormalization point u. Coulomb effects on length scales > 1/fi have 
been removed from it. One would expect the coupling constant Cq(/j.) to be within a few 
percent of the corresponding coupling constants for pn and nn scattering. Our result is 
non-perturbative both in the strong coupling and in the fine structure constant which 
is seen to enter in the combination a In a. This is a consequence of the Coulomb force 
becoming strong when the energy goes to zero. 

Since we are not including pions in the present effective theory, we must have the 
renormalization mass in the range |l/a^,| < fx < m n . Depending on the exact value we 
see that the Coulomb correction can actually become of the same magnitude as the strong 
interaction. In fact, since the scattering length arj, is negative, it can have a very big 
effect on the hadronic scattering length a pp . This is seen from the right-hand side of (|Tl| ) 
which becomes zero when /x ~ m n as shown in Fig. 3. It is not a very physical result, 
corresponding to the hadronic scattering length a pp being infinite. 

The exact numerical properties of our result is obviously dependent on the regular- 
ization scheme we use. As a check, we have also calculated the Coulomb-dressed bubble 
(||) using instead a momentum cutoff A as a regulator. It will then appear in the loga- 
rithmic term of ( pT| ) instead of ll and the numerical constants inside the parenthesis will 
be different. In fact, the new result can be obtained from by letting fi ~ 0.56A. This 
corresponds to a blowup of the hadronic scattering length for A ~ 240 MeV. 



4 



15.0 



10.0 



-5 
I 



5.0 - 



0.0 - 



-5.0 ' — 
20.0 



80.0 



140.0 



200.0 



\X (MeV) 



Figure 3: Dependence of the inverse scattering length on the renormalization point. Dashed curve 
gives result with only Cq interaction, while the solid curve also includes the Ci interaction. 



From a physical point of view such an unstability in the result is not acceptable. 
Instead we should like to see the dependence on the regulator mass to be as small as 
possible. The obvious way to cure this instability is to include higher order interactions 
in the effective Lagrangian ([!]). In the S-wave channel there is only one such interaction, 

C 2 = \c 2 {pa 2 V 2 p)(pa 2P ^ + h.c. (12) 

where the operator V = l/2(^ — V). Kaplan, Savage and Wise[jJ| showed that it can be 
included in a perturbative way and gives directly the effective range ro of the proton. The 
resulting Feynman diagrams can be obtained from the original diagrams by replacing one 
Co vertex in each diagram in Fig. 2 by a C 2 vertex in every possible way. We then find the 
following correction to the scattering amplitude 

when we neglect higher orders in the fine structure constant. It depends crucially on the 
properties of the Coulomb wavefunction at the origin. While 1^(0)1 = C„, the double 
derivative is divergent which is easily seen from the Schrodinger equation. A regulated 
expression can be obtained from the exact wavefunction which we transform to momentum 
space. The resulting integration is then done in the PDS scheme. In this way we find 

i2 I 2 i „ j,, r ,. i 1 „,2 njr2 



r(0)V^(0)\re 9 = ^ ( -j? + aM/i + -a 2 M 2 ) (14) 
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after a rather elaborate calculation. The term proportional to fi is again from a PDS pole. 
Dropping the higher order a 2 term, we thus have the following correction to the inverse 
scattering amplitude 



1 \ e- 2ia °C 2 (n) 



Tsc(p)J " C^C 2 ^) 



p 2 -aMfi) (15) 



Comparing with the effective range expansion (Q) we then obtain for the effective range 
of the proton simply r$ = 8ttC2/MCq. It is therefore not modified by the Coulomb 
interactions. This is consistent with experiments which show that the effective range for 



pp is the same as for pn and nn to within a few percent [12] 



With the coupling constant C2 determined in this way, we can now calculate its 



contribution to the scattering length by taking the limit of zero momentum in (15). The 
PDS pole is seen to give an extra term in the previous result (|ll]) which now is changed 
into 



1 pp(^) a pp 



+ aM 



1 ^ , 1 3 n 1 

ln ^M +1 -2 CE ~r ro 



(16) 



Although this new term has been calculated in perturbation theory, it comes into the 
result with the same general magnitude as the first Coulomb correction and actually 
stabilizes the result when considered as a function of the renormalization point \x. This is 
shown in Fig. 3 where the corresponding curve varies little in the whole interval 80 MeV < 
fi < 200 Mev. The blowup of the scattering length for [i approaching the pion mass 
has now been turned into a rather smooth behaviour in the same region. In fact, we 
find — l/o( / u = m-jr) = 6.6 MeV. For comparison, it should be noted that — l/a pn = 
8.3MeV and -l/a nn = 10.5MeV while (|) gives -l/a pp = 11.6MeV. These differences 
in scattering lengths correspond to much smaller differences in the corresponding coupling 
constants Co- The strong scattering length in our approach is scheme dependent and it 
will not enter explicittly in any physical results. 
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